Probing Dark Matter with AGN Jets 



Mikhail Gorchtein, 1 ^ Stefano Profumo, 2 '[j] and Lorenzo Ubaldi 2 '[*] 

1 Center for the Exploration of Energy and Matter, Indiana University, Bloomington, IN 4^408 
2 Santa Cruz Institute for Particle Physics and Department of Physics, 
University of California, Santa Cruz CA 95064 

We study the possibility of detecting a signature of particle dark matter in the spectrum of gamma- 
ray photons from active galactic nuclei ( AGNs) resulting from the scattering of high-energy particles 
in the AGN jet off of dark matter particles. We consider particle dark matter models in the context 
of both supersymmetry and universal extra-dimensions (UED), and we present the complete lowest- 
order calculation for processes where a photon is emitted in dark matter-electron and/or dark 
matter-proton scattering, where electrons and protons belong to the AGN jet. We find that the 
process is dominated by a resonance whose energy is dictated by the particle spectrum in the 
dark matter sector (neutralino and selectron for the case of supersymmetry, Kaluza-Klein photon 
and electron for UED). The resulting gamma-ray spectrum exhibits a very characteristic spectral 
feature, consisting of a sharp break to a hard power-law behavior. Although the normalization of 
the gamma-ray flux depends strongly on assumptions on both the AGN jet geometry, composition 
and particle spectrum as well as on the particle dark matter model and density distribution, we 
show that for realistic parameters choices, and for two prominent nearby AGNs (Centaurus A and 
M87), the detection of this effect is in principle possible. Finally, we compare our predictions and 
results with recent gamma-ray observations from the Fermi, H.E.S.S. and VERITAS telescopes. 
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I. INTRODUCTION 



The innermost regions of active galactic nuclei (AGNs) correspond to locations where the highest dark matter 
densities in the universe are believed to occur. This remarkable feat, and the fact that AGNs are well-known sources 
of high-energy particles, naturally leads to the question of whether anything could happen and be observed as those 
high-energy particles in the AGN jets traverse high-density dark matter regions. One possibility, first cleverly suggested 
some time ago by Bloom and Wells in Ref. pQ, is that, for AGN jets pointing off of our linc-of-sight, an isotropic photon 
emission might result from the scattering of the high energy electrons in the jet off of dark matter particles. The 
estimates of Ref. [I], unfortunately, and for reasons we will discuss in the present analysis, led to rather pessimistic 
conclusions, and this novel route to search for particle dark matter was basically dismissed. 

With the advent of the Fermi Large Area Telescope (LAT) [2j, and with the amazing recent results from ground- 
based Cherenkov Telescopes, most notably H.E.S.S. [3J, we deem it timely to re-consider the original proposal of 
Bloom and Wells. In particular, we here carry out a more in-depth analysis and we improve on the results of Ref. [T] 
in several respects, including: 

1. a complete calculation of the relevant cross section for neutralinos in the minimal supersymmetric extension of 
the Standard Model and for the lightest Kaluza-Klein particle of Universal Extra-Dimensions 

2. data-driven and semi- analytic models for the dark matter density distribution of the target AGNs 

3. observationally motivated jet models 

4. a novel discussion on proton jets and proton-dark matter scattering 

5. a comparison with the recent Fermi-LAT and H.E.S.S. data. 

Perhaps more crucially than the refinements listed above, though, the main novel theoretical fact we point out here is 
that the electron (or, as we shall also consider here, proton) scattering off of dark matter produces a unique spectral 
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feature, physically associated to a resonance that is present for several motivated particle dark matter models. This 
property of the cross section was not recognized in the seminal work of Ref. [1] due to a simplifying assumption in the 
basically only dimensional estimate of the relevant cross section. The resonant channels we find, quite interestingly, 
dominate the overall size of the scattering cross section, yielding much larger photons fluxes than those estimated in 
Ref. [J. 

Observationally, the predicted spectral feature corresponds to a sharp cut-off in the spectrum of the outgoing 
photons, a sudden spectral break that might in principle be detectable by both Fermi and Cherenkov Telescopes 
observations. The energy corresponding to the cut-off depends non-trivially on the particle spectrum of the dark 
sector, and specifically on two key mass scales, including the mass of the dark matter particle. Whether or not this 
predicted signature will actually manifest itself in gamma-ray observations depends however critically on the details 
of the AGN jets and on the dark matter density distribution, in addition to the dark matter particle model. Also, 
disentangling this feature from off-jet astrophysical gamma radiation in the AGN is a highly non-trivial exercise, 
which we leave to future studies, whereas the focus of the present study is limited to the structure and intensity of 
the exotic dark matter related signal. 

The differential photon flux from electron or proton (e[p]) scattering off of dark matter particles is simply expressed 
as the following integral over the energy of the impinging e[p]'s: 

The integrand is the product of three factors: 

1. the first factor, 5dm, involves the dark matter density profile in the AGN, and is defined in Eq. ^ below; 

2. the second factor depends upon the target AGN, and in particular on the AGN distance cZagn and on the 
spectrum of the electron or proton flux; 

3. the third factor, finally, depends upon the dark matter particle model, and it involves the total cross section 
for e[p]~x scattering (where \ indicates the dark matter particle) into any final state that contains a photon, 
divided by the dark matter particle mass M x (the cross section is computed at a scattering angle 0$ between 
the direction of the AGN jet and the line of sight). 

Conceptually, our analysis addresses, one after the other, the three factors listed above, and eventually leads to 
estimates of the full integral in Eq. ([IJ. For the sake of comparing theoretical estimates with actual gamma-ray 
observations, choices need to be made for both the target AGNs and the particle dark matter models. For the first, 
we focused on two nearby AGNs, M87 and Centaurus A, for which recent Fermi-LAT data were presented in Ref. [H[S]. 
Observations at other wavelengths for these two well-studied objects also allow for educated estimates for the relevant 
parameters for the AGN jets entering the second factor above. As far as dark matter particle models are concerned, 
we rely on two widely studied weakly interacting massive particle models: supersymmetric neutralinos [5] and the 
Kaluza-Klein B^ 1 ' particle from the Universal Extra Dimensional scenario [7J. On a side note, we remind the Reader 
that other dark matter particles, specifically axion-like particles, can also produce observable effects in the gamma-ray 
flux from AGNs, for instance via the attenuation produced by photon to axion-like particle conversion in the presence 
of a magnetic field near the source, as pointed out in the recent studies listed in Ref. [5]. 

The organization of the paper is as follows: In section [n] we discuss the estimate of the dark matter density profile 
and the calculation of the first factor in Eq. ([I]) ; In section III we discuss the spectral distribution of electrons and 



protons in the AGN jet, that enter the second factor of the above integral. We then calculate and study the cross 
section for the process of interest in section [TV] that will allow us to calculate the third and last factor of the integrand 
of Eq. . Finally, we put all three factors together to compute the total photon flux in section [Vj and we discuss 
our results in relation to the sensitivity of current telescopes, and we especially compare our results with the recent 



Fermi-LAT observations [H |S]- We devote sec. VI to a discussion of the case of proton jets, and we conclude in 
sec. |VII| We leave some considerations on relativistic effects in the particle energy distribution in the jet to Appendix 
A and the details of the cross section calculation to Appendix B. 
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II. THE FIRST FACTOR: DARK MATTER DENSITY PROFILE 



The first factor in Eq. ([I]) corresponds to the line-of-sight integral of the dark matter density in the AGN under 
consideration, i.e. to an effective average of the dark matter density times a relevant length-scale associated to the 
AGN jet: 

pro 

5vm =< Pdm-Rdm >= / Pdm(V) dr, (2) 

where the integration runs between r m i n , the minimum distance from the AGN center at which the scattering process 
we study takes place (thus, effectively, the base of the AGN jet), and ro, an upper limit of integration that corresponds 
to the distance at which the AGN jet peters out. Also, in the Equation above, /3dm(?~) is the dark matter density 
profile as a function of the distance from the central black hole. We expect only a mild dependence on the precise 
value of the upper limit of integration r , since dark matter density profiles fall off steeply with radius, whereas it 
will be important for our calculation to model as accurately as possible the innermost portion of the dark matter 
distribution, taking into account the gravitational effects of the central super-massive black hole on the equilibrium 
density distribution. 

The two AGNs we consider in our study, Centaurus A and M87, lie at a estimated distances of 3.7 and 16 Mpc, 
respectively [§1 OH] • Such distances are large enough to make it impossible to have observations that can resolve 
regions as close as 100 or 1000 Schwarzschild radii from the compact object in the AGN core. In addition, dynamical 
tracers of the dark matter density distribution also usually provide very limited information at small radii. We thus 
necessarily need to rely on extrapolations to the dark matter density distribution, which we base on both results of 
N-body simulations and on theoretical studies on the effect on the dark matter density profile of the presence of a 
compact object (the AGN super-massive black hole) and of the surrounding baryons and stars. To this end, we adopt 
the theoretical results presented in Ref. [11] and [12] . 

Working under the assumptions that the dark matter particles are collisionless and that the central black hole 
grew adiabatically by accretion of gas and stars, Gondolo and Silk [TT] found that the dark matter forms a dense 
central spike. The profile they obtain for an initial dark matter density distribution with a power-law inner profile 
PiniM oc t~ 1 reads: 



( \ ( \ P'( r )P 
pDu[r) =p sp [r) = -^y^ 



'core 



(3) 



where 



M x /(Mo%) P'(r) = PR9,{r) f^Y"" ' (4) 



In the Equation above, p R = p 7 , g 1 (r) ~ (l - ^f^) 3 , R Bp = a 7 r (^f) ' \ 7s P = with a 7 a 

numerical coefficient that depends on the initial dark matter density profile slope 7. pdm vanishes for r < 4i?s, 
particles on smaller orbits being accreted on the super-massive black hole [TT]. r = 4i?g thus sets the lowest possible 
value for r min . Some numerical values for quantities that enter the profiles defined above are listed in Table [TJ for 
Centaurus A and M87, together with some nominal assumptions on the particle dark matter models. In particular, in 
Tab. [TJwe consider annihilation rates with an upper bound set by the typical pair-annihilation cross section required 
by obtaining a thermal relic abundance close to the cosmological dark matter density, (av)o ~ 10~ 26 cm _3 s _1 , and 
with a lower bound close to what found in the case of bino-like neutralinos in the MSSM with co-annihilations setting 
the required relic density to the desired value [13], (av)o ~ 10~ 30 cm _3 s _1 . 

A few years after the work of Gondolo and Silk, Gncdin and Primack [12] pointed out that the central regions 
around super-massive black holes usually harbor stars that can gravitationally scatter the dark matter particles and 
cause the distribution to evolve towards an equilibrium solution that differs from what found in Ref. [11] . Taking 
scattering off of stars into account, Gnedin and Primack found a universal inner profile that scales as r~ 3 / 2 , which 
we write as 



3/2 



pDM(r)=po(-) , (5) 



where we set a = 10 5 Rs- Notice that the Gnedin and Primack profile of Eq. ^ is a more conservative choice than 
the one taken in Ref. [TJ, where pdmM ~ r~ 18 . 

The profiles of Equations ^ and ^ are assumed to be valid only in the innermost AGN regions, which are also 
those most relevant to the calculation of the quantity in Eq. ([2]). To normalize the profiles, we need an estimate for 
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the enclosed dark matter mass in the innermost regions around the central point sources. At radii relevant for the 
determination of the black hole masses, typically 10 5 Rs [H1[T5], we require that the enclosed dark matter mass be, 
at most, as large as the uncertainty over the black hole mass, i.e.: 

rW 5 Rs 

/ dr4irr 2 p DM < 3 x 10 7 M Q (6) 

for Centaurus A and 

"10 5 i?. s 

dr47rr 2 p DM < 5 x 10 8 M (7) 

for M87. Note that the lower limit of integration, ri ow = 4i?g is obvious for the profile of Eq. ([3|, while for the profile 
of Eq. ([5]) the lower cutoff will be set by the requirement that the dark matter not annihilate away during the black 
hole lifetime, which corresponds to densities p ~ M x j ' {{(Jv)^tBii)- For the AGN and particle dark matter parameters 
we consider, ri ow — 100i?s is an appropriate choice for the profile by Primack and Gnedin. 

Using the Gondolo and Silk distribution, we notice that one runs in the additional complication that the integrands 
111 ^ and depend non- linearly on po, when determining the latter quantity from Eq. ([6][7| above. We note 



however that the main contribution to the integral comes from r 3> Rs- In that regime we have p colc 3> p'(r) 

so we can approximate pdm (?*) — p'( r )- We can then rewrite p'(r) = p^'' \^ s -) 9~<( r ) ( "F") > where R' sp = 

1 

ot^ro ( M rf" ) 3 '' ■ Now Equations (6) and (7) are linear in po and can thus be easily solved. 

We note that the normalization used here is also consistent to available observational data. As a cross-check, we 
integrated our profiles to calculate the dark matter mass enclosed in 1 kpc and compared it to the mass profiles at 
small radii obtained from astronomical observations and given in Ref. [TO] and Ref. [5] for Centaurus A and M87, 
respectively. In both cases we found the enclosed dark matter mass to be one to two orders of magnitude below the 
total mass at a radius of ~ 1 kpc, thus concluding that our normalization is reasonable and compatible with available 
data. 

We now are in a position to compare the values of the integral <5dm m Eq- ^ - the first factor in the overall flux 
calculation outlined in Eq. (JlJ - corresponding to the two different profiles discussed. We show our results in Fig. 1(a) 



and 1(b) In each one of the two Figures, corresponding to Cen A and to M87, respectively, we show three panels: 
in the upper left panel we plot the dark matter density as a function of the distance from the AGN center. We set 
the original slope (prior to adiabatic compression) of the dark matter density profile 7 = 1. Notice that p'{r) is the 
same for all the curves shown. Also notice how the spike at small r only emerges when the parameters are such that 
Pcorc 3> p'(r), in which case pbm{t) — p'{r). The lower- left panel shows the (very mild) dependence of the 5dm factor 
on 7. Finally, the bottom-right panel illustrates how 5dm depends upon the assumed value of the distance between 
the origin of the AGN jet and the AGN central region, in units of Rs- Again, we set 7=1, which, incidentally, 
corresponds to 7 sp = 7/3. 



The four lines in Fig. 1(a) and 1(b) correspond to different assumptions on the pair-annihilation rate of the dark 
matter particle (which we set to the two representative values (av)o ~ 1CP 26 cm _3 s _1 and {av)o ~ 10 -30 cm _3 s _1 
quoted and motivated above) and on the black hole age £bh (for which we consider typical choices for the ages of 10 8 
yr and of 10 10 yr). The combination ((crv)otBH) sets, for a given dark matter particle mass, the maximal density that 
survives complete annihilation, as can be seen directly from inspection of the upper-left panel. While the four lines 



Centaurus A 



M x [GeV] DM mass 

Mbh [M q ] Black Hole mass 

7?s [pc] Schwarzschild radius 

£bh [yr] Black Hole time 

{<jv)o [cm 3 /s] Annihilation cross section 

Ct-y 

ro [kpc] upper limit of integration 



100 

(5.5 ±3.0) x 10 7 
5 x 10~ 6 



10 8 = 10 1 



10 



-30 



10 



-2(i 



0.1 
15 



M87 
100 

(6.4 ±0.5) x 10 9 
6 x 10~ 4 
10 8 + 10 10 
10" 30 - 10- 26 
0.1 
15 



TABLE I: Values for parameters relevant to the Gondolo & Silk dark matter profile. The black hole mass is from Ref. P3] for 
Centaurus A, and from Ref. [15] for M87. 
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FIG. 1: Top. The dark matter density profile near the central AGN regions, for 7=1 and for four different choices for the 
parameters < av >o and tsH- Bottom. We show how the integral Sum depends on the slope of the density profile prior to 
adiabatic compression 7 (left panel), and on the lower limit of integration r m i n , in units of Rs (right panel). Figure (a) refers 
to Centaurus A, while (b) to M87. 
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all correspond to the density profiles described by Gondolo and Silk, where the density of dark matter very close to 
the core is dominated by p CO ic, resorting to the profile by Gnedin and Primack we get values very close to the lowest 
curves showed in Figures 1(a) and 1(b) Notice that a profile such as that employed in Ref. [1] would have fallen in 
between the lines we show in those figures. 

We notice that the most promising dark matter models are those with a low pair-annihilation cross section: this is a 
quite remarkable feat, as this is the exact opposite of what is found for any other indirect dark matter detection method! 
Here, the reason is that the less the dark matter pair-annihilates, the more particles can be densely concentrated around 
the central regions of the AGN, and the more probable the scattering of electrons and protons off of these particles is. 
Also, the younger the central black hole, the better, since fewer particles have annihilated away since the black hole 
birth. Interestingly, the dependence on the distance of the base of the jet from the center of the AGN is relatively 
mild: only for extremely spiky profiles, such as the one represented by the blue line, do we find an order-of-magnitudc 
dependence on r m i n , if the latter is within a few hundreds Rs, as is usually expected. Also, the slope of the dark 
matter density profile before adiabatic compression onto the black hole is virtually irrelevant, up to a factor 2 (lower- 
left panel) . We also point out that the final value of the relevant parameter <5dm is mainly dictated by the innermost 
regions of the dark matter density profile. 



Comparing Cen A (Fig. 1(a)) with M87 (1(b)), we find that for the latter, featuring a significantly more massive 
central super-massive black hole, and hence a much larger Rs, the dependence on r m ; n is somewhat stiffer. Also, as 
apparent from the upper-left panel, in some cases the spike is cut off by the r > 4i?s limit beyond which particles 
infall in the black hole [TT]. Overall, for the same particle physics model, a reasonable range of expected values for 
5dm for Centaurus A falls between 10 8 and 10 11 in units of solar masses per pc 2 , while a similar but slightly narrower 
range is expected for M87 (3 x 10 8 < S DM /(M Q /pc 2 ) < 10 11 ). 



III. THE SECOND FACTOR: HIGH ENERGY PARTICLES IN THE JET 

The second factor in Eq. ([!]) depends upon the nature of the high-energy particle content of the AGN jet, including 
its composition, geometry, spectrum and overall energy content ("luminosity"). We assume that the jet consists of 
either or both relativistic electrons and protons. For our purposes, the jet geometry is largely irrelevant: as we noted 
above, the exact position of the base of the jet is not fundamental for our purposes as long as it lies within ~ 100 Rs, 
which is believed to be the case for the AGNs under consideration here [H]. The reason is that, as we saw in the 
previous section, the Sum integral is not much affected by a variation of the lower limit of integration r m i n , as long 
as r min < 100 Rs- Other details of the jet geometry, such as the opening angle, for example, are also not crucial to 
our calculation. One exception is the inclination of the jet with respect to the line of sight, which we will discuss in 
the next section, and which enters the angle at which the cross section is calculated. The crucial information we are 
after is, rather, the energy distributions of electrons and protons that enter the second factor in Eq. (Til. 

Let us focus on the case of electrons first. We adopt here the commonly used "blob geometry" |17j: the electrons 
move isotropically in the blob frame with a power law energy distribution, and the blob itself moves with respect 
to the central black hole with a moderate bulk Lorentz factor Tb = (1 — /^fj) -1 ^ 2 - For Centaurus A and M87, Tb 
is estimated to be ~ 3, although this depends on assumptions on e.g. the detailed gamma-ray emission mechanism 
(see e.g. the discussion in Ref. [5])- In this section, primed quantities refer to the blob rest frame, while unprimed 
quantities refer to the black hole rest frame. Motivated by jet models employed to reproduce the results of gamma-ray 
observations [U [5j, we use a broken power law energy distribution for the relativistic particles in the jet: 



d$ AGN 1 



2_ 

7br 



for 7min < l' < 7max, (8) 



where 7' = E' /m e , and the constant k e can be evaluated from the kinetic power of the jet, L e . First we need to boost 
this distribution to the black hole frame (details of this calculation are presented in Appendix [A]) . We find 

1 fTm** W<J>AGN pi j f1~x d$ AGN 

<W <*/^b-(V)= / twt^V^ / *y^f— ( 7 r s (i-W), (9) 



where \i! = cos 8' , fi = cos 9, 9' is the polar angle of the electrons in the blob frame, 9 is the polar angle in the black 
hole frame. The distribution on the left hand side is isotropic, as we mentioned, and does not depend on //, but the 
boosted distribution, on the right hand side, depends non trivially on /i. Now we can evaluate the constant k e using 



Le= r li r W m rf "- f— ( 7 r B (i- M). (10) 
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The values of the parameters that appear in the Equations above are determined observationally and are taken here 
from the analyses of Ref. [5] for Cen A and of Ref. [4] for M87. They are summarized in Table [TT] 

Notice that for the case of Cen A, in Ref. [5] Abdo et al. point out how it is not possible to fit the Fermi and the 
H.E.S.S. data simultaneously using a single-zone synchrotron/synchrotron self Compton (SSC) model. One possibility 
mentioned there is that the H.E.S.S. emission originates from a different part of the jet than that responsible for the 
gamma rays observed by LAT. Our study suggests that another possibility, electron-dark matter scattering, is in 
principle viable to explain at least part of the Fermi gamma-ray detection. Thus, when choosing the appropriate 
normalization, we are not strictly bound to using a set of parameters that gives the best fit under the assumption 
of SSC. Looking, as a guidance, at the four different sets of parameters provided in Ref. [5], we employ the one that 
gives the best fit to the H.E.S.S. observations (since, with hindsight, the latter cannot be fitted by the emission model 
we consider here), but we use a slightly higher value for the jet power. As stated in Ref. [5], the jet powers quoted 
there are very conservative and can be considered as lower limits: It is thus safe to explore even higher values (and 
we will later on) as long as they are below the Eddington luminosity which, for CenA, is ~ 10 46 erg/s. 

Interestingly, we find that for ~ 3 most of the electrons are seen to move in a narrow cone around the direction 
of the jet in the black hole frame (see Appendix |A| for further details). When we compute the photon flux from 
Eq. 0), 

we will use for the second factor 

' ' ( 7 r(i-^)). (ii) 



rf AGN dE e ^AGN JQ.9 F s(l " PbP) dj 

Restricting the region of integration over /.i from 0.9 to 1, one loses roughly 20% of the electrons, while this choice 
restricts the high-energy scatterers to a quite collimated jet. The angle of the emitted photon with respect to the 
jet can now be considered to have an uncertainty of around ± arccos(0.9), comparable to (or much smaller than) the 
astrophysical uncertainties on the angle that the jet has with the line of sight, and therefore not a major concern for 
the accuracy of our final conclusions. 



IV. THE THIRD FACTOR: CROSS SECTION 



In this section we compute the cross section for the process "electron + dark matter — > electron + dark matter + 
photon" using two different dark matter particle models: (1) the Minimal Supersymmetric Standard Model (MSSM), 
where the dark matter candidate is assumed to be the lightest neutralino, x (2) Universal Extra Dimensions 
(UED), where the dark matter is the lightest Kaluza-Klein particle (LKP), assumed to be the B^ (first Kaluza-Klein 
mode of the hypercharge gauge boson B) [7]. The structure of the calculation is as follows: We first perform a 
full calculation for the supersymmetric model. We then infer from that calculation which terms give the dominant 
contributions to the cross section, and we apply that knowledge to the second model, UED, for which we only compute 
the analogous important terms. 



A. Neutralino Dark Matter 



We study the scattering process of electrons off of neutralinos with the emission of photons at the lowest significant 
order, i.e., symbolically, e + x^ e + X + 7- I n principle there are several contributing Feynman diagrams, but without 



Parameter Symbol Cen A M87 

Low-Energy Electron spectral index si 1.8 1.6 

High-Energy Electron spectral index S2 3.5 3.6 

Minimum electron Lorentz factor 7min 8 x 10 2 8 x 10 2 

Break electron Lorentz factor 7b r 4 x 10 5 4 x 10 3 

Maximum electron Lorentz factor 7max 10 8 10 7 

Jet power in electrons [erg s" 1 ] L e 3 x 10 43 10 44 



TABLE II: Values for parameters in the electron energy distribution. For M87 the values are taken from Ref. [4, while for 
CenA they are from Ref. 5 . Notice that in Ref. [5] four different sets of parameters are used for four different fits. The one 
set we show here is the one that produces a good fit to the H.E.S.S. data, but, as explained in the text, we use a higher value 
for the jet power. |3j.) 



e(k> 



x(p)- 



J 7(?) 




e(fc') e(Jfc) 



x(p') x(p) 




e(ifc') e(fc) 



x(p') x(p) 





(a)Mj (b)M 2 (c)M 3 

FIG. 2: Feynman diagrams for the process e(k) + x(p) ~^ e {k') + xip') + "/(l)- 



loss of generality and accuracy we can restrict our attention to the s-channel diagrams shown in Fig. [2] In fact, when 
the exchanged scalar electron (selectron) is on shell we have a resonance and these three diagrams are dominant, so 
it is fair to neglect the diagrams in the t- or u- channels, which never feature an on-shell particle in the intermediate 
state. Notice that no resonances at all appear from the two amplitudes considered in Ref. pQ - one of the crucial 
differences between the present results and the findings of that study. By only considering the process ex — > §7, in 
fact, the intermediate state selectron can never, kinematically, be on-shell, and thus the process is never resonant. 
Referring to the momenta as labeled in Fig. [2j we adopt the following notation: 



t 
t' 

s 
s' 

n, 
a,. 



(k-q) 2 
(k' + qf 
(p + k) 2 
(p' + k'f 
1 



s - Mf - i^Y 
1 



(12) 



2 \ (Aff-M 2 ) 2 



is the selectron energy-independent width 2 , with M x 



where Mg is the selectron mass 1 and T = (a 2 L + aj^) 
the mass of the neutralino. 

We present a detailed calculation of the squared amplitude | J2\ spins l-^| 2 m AppendixlBl In this section, instead, 
we focus on those terms that, due to special conditions, feature large enhancements. Theseterms correspond to: 

(i) the resonance for s/s — Mg, i.e. when the selectron is on-shell; 

(ii) a logarithmic enhancement, when the photon is co-linear with the final electron. 
The terms in the cross sections that feature both (i) and (ii) are (see Appendix Ib| : 



A, spins 



KWM 



AH 2 _ 
2 I — 



2e 2 (ai + a i J^f{pk){p' q ), 



\ {e;(\)M^ u (X)Mr+h.c} = - 2e 2 (4+ 4)^111,, | 2 W 



Mi 



'F2 



sT 



A, spins 



}(p'k')4(pk) 2 . (13) 



Note that the enhancement (i) comes from |II S | 2 , when s = Mf , while (ii) comes from 1/t'. Besides (i), there is another 
resonance, associated with n s / for s/s' = Mg, which depends on the final photon energy. This latter resonance is 
always less important than the first one (see e.g. Fig. [3]). Notice that one might expect that we should include 
the term ZRe[M.2-M.^\, which is also logarithmically enhanced. However that term features a factor Re^^n*,) which 
produces a much smaller enhancement than |n s | 2 . Thus, the 2Re[A^2A / li] term is less important than the ones 
considered in Eq. (13) and we can safely neglect it. Our numerical work supports this claim, as we will also see e.g. 
in Fig. [5] below. 



1 For simplicity, we assume the right-handed and the left-handed selectrons to be degenerate in mass. 

2 See the appendix for the definition of the couplings and a^j. 
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We can now compute the overall differential cross section 



da 



{2Ti) i 8 i {p + k-p' -k' - q)- 



1 



d 3 k' 



d 3 p' 



4 ^ 

X, spins 



AM X E 2E'{2n) 3 2£ 7 (2tt) 3 2E' n (2tt) 3 
{\e;(\)M%\ 2 + {e;(X)M^e u (\)Mr + h.c}) . 



(14) 



We first concern ourselves with calculating the kinematics. We work in the approximation of massless electrons, 
except when we calculate logarithmic terms (see Eq. (17l). We assume the neutralino to be at rest, p = (M x , 0,0,0), 
and we will eventually integrate over the spectrum of incoming electrons. The quantity that is actually measured, the 
photon energy, in our notation is q. Therefore, we need to integrate over the p' , k! phase space. The p' integration is 
done with the ^-function, so only k' remains. The scalar products in Eq. (131, 



{p'k') = (p + k-k' -q,k')-> (pk) - {p + k, q) 
(p'q) = (p + k - k' - q,q) -)• (p + k,q) 



(15) 



do not depend on angles involving k' . In order to have the logarithmic enhancement, we want the photon to be co- 
linear with the final electron: this corresponds to the so-called co-linear approximation, that will also further simplify 
the calculation of the kinematics. After integrating over d 3 p' with the delta function, we are left with the energy delta 
function that gives us the on-shell condition for the final electron energy: 6(En + E — E 1 — £" — E' N ) = S(E' ' + E 1 — E, 
with the usual Compton laboratory-frame energy Eq 

and final photon. Finally, 



^c) 



E 



l+#-(l-cos( 



6 being the angle between the initial electron 



1 



1 



dE 7 dQ 



{ 2,f,2El 2e ^ + a ^ 
Er 1 {M x + E{l-cosO))- 

dflk' ■ 



in., 



s' - Ml 



v^T 2 



4EM X {EM X - E 1 (M x + E(l ~ cos 0))) 



f 



(16) 



and the integral over the solid angle can easily be evaluated (here one has to keep the electron mass finite) to give 

'4E' 2 



7r In 



(17) 



Above, E' M = Af v + r^^ai is the energy of the neutralino in the final state. Note that from the delta function 

5 jv X M x +E(l— cos 0) GJ 



we have the condition E' 



E i - TT~- 
x i i 



E 



- (1 — cos i 



which leads to a lower limit on the initial electron energy for a given 



angle between the jet and the observer, cos0, and outgoing photon energy E 7 



E,, 



E. 



7 



1 - |^(1 - cos 6>) 

A few comments are in due order at this point: 
• For a given photon kinematics, only those particles in the initial electron spectrum with E > E„ 



(18) 



contribute. 



From Eq. ( 18 ) we see that the photon energy cannot exceed ^ 



AL 



• The cross section in Eq. ( 16 1 contains two resonances: one for s = Mf , that doesn't depend on _E 7 , and a second 
one for s' = M~, that instead does depend on E 1 . The latter, however, gives a peak which is always significantly 
lower than the former. 

• When we scan over photon energies, E m i n gets larger as E 1 increases. Eventually E m i n gets big enough that 
s > M~ and we lose the first resonance. This happens when 



E-, 



M X (M§ 



Ml) 



2MI + (Mf - Af2)(l - cos 9) ' 



(19) 
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FIG. 3: The differential cross section of Eq. (BIO I as a function of the initial electron energy E for different values of the final 



photon energy i? 7 . Note that we use for this Figure the full result as calculated in Appendix, to which Eq. (161 is a good 
approximation . 
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FIG. 4: Feynman diagrams for the process e(k) + £> (1) (p) — > e(fc') + £> (1) (p') + 7(g). 



Some of these features are illustrated in the plots of Fig. |3j where we show the full differential cross section of Eq. ( B10 ) 
as a function of the initial electron energy E for different values of the final photon energy E 1 (not the approximated 
expression we give in Eq. (16)). In the first plot, where we set E 1 = 10 GeV we see the two resonances corresponding 
to s = Mf and s' = Mf . In the second plot, where E 1 — 30 GeV, we still have both resonances, but the line starts 
at a larger value of E = E m i n because we of the increased value of E y . Finally, in the third plot, with E 1 = 50 GeV, 
Ernin is above the value needed for the first resonance to occur, and we only see the second one. As explained above, 
the first resonance does not move when we vary E 1 , while the second one does. Note also that, when present, the 
first peak is much higher than the second one. For all three plots we fixed the following parameters: Mg = 100 GeV, 
M x = 60 GeV, 9 = 68° (the latter angle corresponding to the line-of-sight angle of the Cen A jet). 

The Reader may notice the apparent difference between the present calculation and the one presented in the work 



of Bloom and Wells, Ref. PQ: while the result of Ref. [T] is of order a 2 L R , ours is 



l l.r- 



However, this is only 



seemingly a discrepancy because the two extra powers of cil,r are cancelled by the selectron width in the denominator 
(e — » e + x is the only decay channel) - see e.g. Eq. (21) below. This further leads to the resonant enhancement 

\n(E 2 /m%), is absent in Bloom and Wells [1. 



l/(mg — m x)- This mechanism, along with the log enhancement 



B. Lightest KK particle Dark Matter 

With the insight gained from the calculation of the full differential cross section for the supersymmetric case, we 
can now turn to the case of UED, where the dark matter candidate is a massive vector boson corresponding to the 
Kaluza-Klein first excitation of the hypercharge gauge boson, a particle that is also assumed to be the Lightest 

KK Particle (LKP) [7]. The relevant Feynman diagrams for this model, where the scalar-electron is replaced by the 
Kaluza-Klein electron, are shown in Fig. |4j 
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We define 



M = M Bm 
8 = M e( i) - M Bm « M. 



(20) 



Note that we assume the KK first modes to be nearly degenerate [7j. Under this assumption, we have a very small 

width 3 r = -^jggl (Y£ + Y%) for the decay of the e (1) , so we are looking for a <5-like peak in our calculation, and we 
can perform the following substitution: 



I 2M e)1) 2 

l P 2 -M 2 -iM eW T ] 



1 

{E-5) 2 +T 2 /A 



kT 



6{E-S). 



(21) 



To compute the differential cross section, we can apply here the same line of reasoning as for the supersymmctric 
case. Referring to the diagrams of Fig. [Z| we only need to compute ~ J2\. sp ins 1-^2 1 2 and | J2x, spm s {M 2 Ml + h.c} 
for the leading contribution. We have 



ieg 2 u(k')j° 



4 „YZW+lf)+Y L Y R M ell) 



iegf 



2t'(E-5-iT/2) 



(k' + q) 2 -ml' (p + k) 2 - M 2 (1) - iMT 

, u(k')(2k' a + 7 < W[*At'+ Y L Y R }^P A u(k), 



M% va 



iegiu{k'){Y R P L + YlPrW — 



+ ¥' + M e(1) 



Tf- 



M. 



e<l) 



-^{Y R P R + Y L P L )u(k) 



ieg\ 



iT/2){E -5 - iT/2) 



(p' + k') 2 - M 2 (l) - iMT ' (p + k) 2 - M 2 (1) - iMT 

u{k'){Y R P L + Y L P R )YW+ l)TW+ l)r(Y R P R + Y L P L )u(k), 



(22) 



with A = R, L the chirality of the electron, and = = (1, 0) in the "laboratory" frame. The resulting spin- 
averaged squared amplitudes are: 



A, spins 



1 + SE cosfl 



9(Y 2 + Y 2 ) 2 

1 9Rn[A/f A/f*l - 1r 2 ^(V 2 I v 2 ) 2 {E — E 1 — &)EE' (E' + E^) 

- ^ 2Re[M 2 M 3 \ - -3e 9l (Y L +Y R ) — — _ 



A, spins 



t'[(E - S) 2 + T 2 /4}[(E -E y - S) 2 + T 2 /4] 



HYj-Y 2 ) 2 



cos 9 



(23) 



Note that we have two scales in this calculation, E ~ E' ~ E 1 ~ <5 <C M. To derive the Equations above we dropped 
terms proportional to S/M. The kinematics is the same as in the neutralino case. Because of the hierarchy between 
the two scales, the calculation is here even simpler. For instance, the Compton laboratory energy is Eq — E, the 
energy of the final B^> is E' B ~ M and the lower limit onto the initial electron energy just reduces to E m i n ~ Ej. 
Finally, the differential cross section for the UED case reads: 



d 2 c 



e 2 9t 



3tt {Y 2 + Y 2 ) 2 + l{Y 2 -Y 2 ) 2 cost 



dE 7 dn q (2tt) 5 UME' b 



{E-5) 2 + T 2 /A 



(E-E,-S)E'(E' + E 1 ) 
1 {E-E 7 -5) 2 +T 2 /A 



In 



AET 
ml 



(24) 



V. PHOTON FLUX 

We now have all the ingredients to compute the differential flux of Eq. ([!]), with educated guesses for the two 
AGNs Centaurus A and M87, and we can thus proceed to compare the process under investigation here with actual 
gamma-ray observations. Prior to comparing with data, we carry out a parallel study of the supersymmetric versus 



3 Fl(Yr) is the hypercharge of the left-(right-)handed electron, g\ is the hypercharge gauge coupling. 
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FIG. 5: We plot Alat-E^^ 1 vs E 1 for the two scenarios we studied: neutralino dark matter (top) and LKP dark matter 



(bottom). To obtain the black line in the top plot we used the full result Eq. (BIO I for the differential cross section. Note how 
the blue line coincides with the black one until the curves drop down. This confirms the validity of the approximation we made 
in the calculation of the differential cross section, namely that the main contribution comes from \M 2 \ 2 and 2Re[X 2 Xg]. As 
expected, the curves drop down when we lose the main resonance in the cross section (see Eq. (19 1 for the SUSY case). The 
two cases give an almost identical signal, apart from a constant that can be traced back to the different coupling constants. 
In the top plot we also show with a yellow line the signal one would obtain considering only the process with no resonances 
studied originally in Ref. pQ. The angle 9 is here fixed at 68°, a value appropriate for Cen A. 



UED dark matter paradigms as far as e[p]-dark matter scattering is concerned. In Fig. [5] we plot the results for the 
supersymmctric case (top) and the UED case (bottom), for comparable values of the mass parameters (dark matter 
particle mass set to around 300 GeV, and mass splitting between the dark matter particle and the intermediate state 
scalar-electron or Kaluza-Klein electron of 5 GeV; notice that in the case of supersymmetry while such a small mass 
splitting might not be natural, it is sometimes warranted by the requirement of a co-annihilation amplitude to enforce 
the correct thermal relic abundance for the lightest neutralino) . The yellow curve shows the full result for the ex — > §7 
amplitudes considered in Ref. pQ (not the approximated estimate to the cross section given there) i.e. with the final 
selectron on shell, and illustrates one of the reasons why the final fluxes calculated here are significantly larger than 
those presented in that study. 

Note that the curves for both UED and supersymmetry are very similar. The fact that the UED curve is slightly 
larger is accidental, and due to the different values of the couplings at the vertices of the Feynman diagrams and 
not to the intrinsic structure of the relevant scattering process. From the plot at the top we also have an explicit 
confirmation of the line of reasoning we outlined in the previous section: the dominant contribution to the process of 
interest here stems from |A^2| 2 and 2Re[A / J2-^3]. The other contributions are not important until photon energies 
become greater than 5 GeV. As we can see, however, the signal drops anyway in that region, so those contributions 
are never crucial. 

We are now in the position to answer the question of whether current experiments might detect a signal coming 
from the process under present investigation. We first consider the case of Cen A [5], and we start focusing on the 
supersymmetric case. In Fig. 6 we plot the spectral energy distribution, i.e. x versus the photon energy, and 

we compare with Fermi-LAT^5] and H.E.S.S. j3] observations. As explained in previous sections, it is reasonable to 
consider a fairly generous range of plausible values for the 5dm integral as well as for the luminosity of the jet, due 
to uncertainties in the physics of the AGN jet and on the dark matter density distribution in the innermost AGN 
regions. In Fig. [6] we show lines corresponding to two different choices for the jet luminosity, that given the results of 
[5] can be regarded as relatively conservative, and we use the most conservative value for 5dm, corresponding to the 



black line in Fig. |l(a) 
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FIG. 6: We show the spectral energy distribution vS v (equivalent to x -j^) versus the photon energy for CenA, for 
four different choices of the neutralino and lightest selectron mass in a supersymmetric scenario. As explained in the text 
astrophysical uncertainties affect both the value of the integral 5dm and the jet luminosity. Fermi-LAT data are from Ref. [5]. 
The H.E.S.S. data are from Ref. [3J. 



The particle dark matter mass is set to 60, 95, 300 and 500 GeV, from top-left to bottom-right, with a range of 
mass splittings between the selectrons and neutralino. Intriguingly enough, we find that the resulting gamma-ray flux 
is detectable for small enough values of the mass (less than about 300 GeV) and a luminous enough jet, even for the 
most conservative assumption on the dark matter density distribution. The spectral structure of the predicted flux of 
photons from electron-dark matter scattering is very hard and has a dramatic drop-off at the edge of the resonance. 
The location of such a sharp break is between 5 and 50 GeV for the models under consideration, the ideal energy range 
for Fermi-LAT observations. We also find that no photons are predicted to be observable by Atmospheric Cherenkov 
Telescopes (ACTs) such as H.E.S.S. It is somewhat fascinating that the Fermi-LAT data from the Cen A system do 
show a hardening of the spectrum right before an energy of 10 GeV. 

Models where the mass splitting between the selectron and the neutralino is much larger than the ones considered, 
are not very interesting here. In those cases, the drop-off of the photon flux happens at higher energy, and one might 
hope to get in the region probed by H.E.S.S. The signal would be too low to be detectable though, mostly due to the 
fact that, in order to hit the resonance, we would now need electrons with much higher energies, that would put us 
in the tail of the electron energy distribution. 

Given the very large values of dark matter density, a potential additional source of gamma rays is direct dark matter 
annihilation close to the AGN center. The flux of such photons can be estimated as 



'7 \ uu l" 



<h\> \ d.V. (av) 

dEl) annihilation _ ^ 2M^d\ & , 



dr 47rrX M (r), (25) 



where the factor dN 1 /AE 1 indicates the differential photon spectrum per annihilation event for the particular dark 
matter model under consideration. Notice that the flux above (and in particular the spectrum encoded in d7V 7 /di? 7 ) 
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FIG. 7: We show the spectral energy distribution vS v (equivalent to E 7 x jg 1 ) versus the photon energy for CenA, for two 
different choices of the masses in a UED scenario. Fermi data are from Ref. [5], while H.E.S.S. data are from Ref. [3]. 



is entirely disconnected (in a model independent way) from the flux originating from e[p]-dark matter scattering we 
discuss here. In principle, one can have the latter without having any flux from the former. Indeed, if we consider a 

small enough annihilation cross section, (crw)o ~ 1CP 30 cm 3 s _1 , typical values of I 4g- j will be ~ 10~ 13 , 

1 / annihilation 

10~ n erg s~ 4 cm -2 . These are comparable to the photon fluxes shown in Fig. |6j coming from the scattering of 
electrons off the neutralinos. Notice that the signal shown in the Figure is peaked at E~ ~ a few GeV, for most of 
the cases considered, whereas the signal from dark matter annihilation would be peaked at Ej < M x , so, even if the 
values of the fluxes were similar, we would see two distinct signatures. Also, in scenarios where the neutralino and 
the selectron are nearly degenerate, the coannihilations of dark matter with selectrons will favor a small value for the 
pair annihilation cross section. This, in turn, would enhance the value of the <$dm and, therefore, the flux of photons 
from the scattering process we studied, while it would result in a less important flux for photons coming from the 
pair annihilation. 

In Fig. [7] we show the results for the UED case, again for the Cen A system, and for the same choices of AGN 
parameters as for the previous Figure. This time the dark matter mass was set to 300 (left) and 500 GeV (right), 
with a mass splitting of 5 (left) and 10 GeV (right), close to what expected in the context of UED [T. Thanks to a 
larger scattering cross section, the UED case gives a potentially detectable signal, for a large enough jet luminosity, 
even with a dark matter mass of 300 GeV. 

The next Figure [8] shows how the Fermi-LAT data from Cen A in the few GeV range can be in principle fitted with 
the emission originating from electron-dark matter scattering, for selected values of the particle dark matter mass, 
and of the astrophysical parameters 5dm and the jet luminosity L, for both supersymmetry (left) and UED (right). 
In both panels, we set the Cen A jet luminosity to the somewhat large value of 3 x 10 44 erg/s. This value can be 
part ly tra ded off for a larger value of <5dm, which we fix here to lie towards the lower end of the expected range (see 



Fig. 1(a)) to 5dm = 2 x 10 s (left) and 6 x 10 s (right) M /pc . The predicted spectral drop-off due to the location 



of the resonance (set by the splitting between the dark matter particle and the intermediate charged particle) follows 
the Fermi data for the two particular choices we made for the dark sector spectrum: 95 GeV for the neutralino mass 
and 110 GeV for the selectron in the case of supersymmetry, and 300 and 315 GeV for the Kaluza-Klein "photon" 
and electron for UED. 

We move on to the case of the M87 jet in Fig. [9j where we show the results for the SUSY case, for identical values 
of the neutralino and selectron masses as in the analogue Cen A case shown in Fig. [6j We compare the predictions 
for the electron-dark matter scattering discussed here with data from Fermi-LAT [J as well as from H.E.S.S [TS] and 
VERITAS [H] (we do not show here comparable data from the Magic Telescope, see [ID])- The H.E.S.S. measurements 
we show are relative to a low state in 2004 (red) and to a high state in 2005 (magenta). In the case of M87, we fix <5dm 
to 10 8 M Q pc -2 , and we show two values for the jet luminosity, respectively 10 44 and 10 45 erg s _1 . Only for a very 
favorable dark matter spectrum, with the dark matter particle featuring a mass of 95 GeV and the selectron mass 
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FIG. 8: We illustrate in this Figure that, with astrophysical parameters adjusted ad hoc, one can fit the data from Fermi-LAT 
for appropriate masses within both a supersymmetric (left) and universal-extra-dimensional (right) dark matter model. 



at 100 GeV, do we find that the process we consider here (again, for the assumed values of the dark matter density 
profile and jet luminosity) is relevant with respect to the observed gamma-ray intensity in the Fermi energy range. 
Also in analogy with the case of Cen A, the spectral drop-off occurs at energies much smaller than those probed by 
ACT's. 

Finally, in Fig. 10 we show the results for the UED case, for M87, with the same jet and dark matter density profile 
parameters as for the previous Figure, and with the same particle dark matter setup as in Fig. The elect ron-LKP 
scattering process does give rise here to a detectable signal, only slightly smaller than the actual Fermi measurement 
at an energy of 10 GeV, for an LKP mass of 300 GeV and a mass splitting with the KK electron of 5 GeV. Increasing 
the mass to 500 GeV reduces the signal significantly (right panel). 

In conclusion, we find that for reasonable assumptions on the AGN jet parameters and on the dark matter density 
distribution around the central massive object of the Centaurus A and M87 AGNs, the Fermi Telescope is in principle 
sensitive to the detection of photons resulting from the scattering of AGN jet electrons off of dark matter particles. 
This conclusion holds for specific choices of the dark sector particle spectrum, but is general enough to include both 
the case of supersymmetric and universal extra dimensional dark matter. 



A. Gamma Ray absorption 

Electron-positron pair production in the process 77 — > e + e~ limits the minimum radius for a compact source 
below which the source is effectively obscured: gamma rays are lost to electron-positron pair production as they 
traverse photon fields surrounding the compact source. Quantitatively, this occurs if r 77 ~ n 7 er 77 i? > 1, where n 7 is 
the target photon number density, cr 77 is the relevant pair-production cross section, and R is the distance from the 
source. Indicating with L 7 the source luminosity at the energy E corresponding to the peak of the 77 — > e + e~ cross 
section for a given gamma-ray energy, we have that this condition translates into the following requirement on the 
photon density: 

n 7 a L l , (26) 

where the numerical coefficient in front of the Equation above depends on the emission geometry. For the gamma ray 
energies of interest to us, E lies in the range between 0.1 and 10 keV. Ref. [21] estimates that the total luminosity from 
the nuclear region of Cen A is ~ 10 43 ergs/s. This implies that the source is transparent to gamma rays in the energy 
range of interest to us, since the radius R at which r 77 ~ 1 is much smaller than R$- In the case of M87, instead, 
Ref. [22] finds that integrating the continuum spectrum of the nuclear region of M87 between 100 /im (~ 1.2 x 10~ 5 
keV) and 10 keV the total luminosity in the highest state is on the order of 3 x 10 42 ergs/s. Since Rs is much larger 
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FIG. 9: The spectral energy distribution versus the photon energy for M87, for four different choices of the masses in a SUSY 
scenario. The Fermi data are from Ref. [2]. The H.E.S.S. measurements [T8] are relative to a low state in 2004 (red) and to a 
high state in 2005 (magenta). The VERITAS measurements are from Ref. [19] . 



for M87 than for Cen A, this implies that also for M87 the in-situ absorption of gamma rays is negligible for the radii 
of interest for the present study. 



VI. PROTON JETS 



So far we have only considered electrons in the jet. Protons are most likely present in the jet as well, and in some 
models that have gained some attention recently, they can even be the dominant component [IT] . In this section we 
analyze the scattering of protons off of dark matter, with the emission of photons in the final state. Because these 
photons are quite hard (> 1 GeV), they will be emitted in the final state by the quark before hadronization (notice 
that they can also be emitted from the squark exchanged in the s-channel). The most straightforward approach 
to compute this process is to study it at the parton level, neglecting details of the hadronization process. We will 
then have the same Feynman diagrams that we already computed in previous sections, with the electrons replaced 
by quarks, and the selectron replaced by a squark. If the incoming proton has momentum fc, the quark will carry 
momentum xk, with < x < 1. The probability of having a quark carrying a given fraction of the proton's momentum 
is set by the proton's Parton Distribution Functions (PDF). The proton PDF indicate that, in general, gluons can 
carry a very large fraction of the proton momentum. One might then worry that we should also consider Feynman 
diagrams with a gluon in the initial state (see e.g. Fig. [TTj). All such diagrams, however, can only be in the i-channel, 
so they can never undergo a resonant enhancement. Therefore, their contribution to the cross section is negligible 
compared to the s-channel diagrams that have resonant behaviors, and we can neglect them for our current purposes. 

The differential cross section, in the co-linear approximation, is then obtained from the one already computed for 
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FIG. 10: We show the spectral energy distribution vS v (equivalent to E% x < ^ L ) versus the photon energy for M87, for two 
different choices of the masses in a UED scenario. As explained in the text, due to astrophysical uncertainties, one can play 
with the value of the integral <5dm an d the jet luminosity within reasonable ranges. The Fermi data are from Ref. |4|. The 
H.E.S.S. measurements [TB] are relative to a low state in 2004 (red) and to a high state in 2005 (magenta). The VERITAS 
measurements are from Ref. [19] . 




FIG. 11: A possible diagram with a gluon in the initial state, 
like this in our calculation. 



As explained in the text, we do not include f-channel diagrams 



electrons by replacing k (the momentum of the electron) with xk (the momentum of the relevant quark, expressed as 
a fraction of the proton momentum) , the mass of the selectron with that of the squark, the mass of the electron in 
the log with a cutoff of the order the QCD scale (we use here 100 MeV). The flux of photons is then given by 



dE 7 J J Xmin ^ d K> \M x dQdE, ) cose \d\ GN dE ) 



'DM j 



(27) 



where fi(x) are the relevant PDF. (For our numerical calculations we used the Mathematica package provided by the 
CTEQ collaboration [23].) 

As far as the AGN proton jet is concerned, we consider here two theoretically and observationally motivated distinct 
cases [T7] : 

1. A jet consists of mono-energetic protons, which is believed to be the case for the innermost AGN jet regions 
[I?] . If the proton energy is E p we thus have 

d$ AGN 

—E— = kl 6(E-E p ). (28) 

The expression above has to be normalized using the kinetic power of the protons in the jet, as deduced from 
observations 

f d$AGN 

L P = dEE—^, (29) 
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which simply gives k\ — L p /E p . 

2. The protons in the AGN jet feature a power law energy distribution. At some distance from the central compact 
source, the protons undergo a shock, which, assuming a simple Fermi acceleration mechanism, results in a protons 
spectrum proportional to E~ 2 . In this case, one has 



p — u I p 



dE V E 

Using the kinetic power to normalize the distribution, we finally get 



(30) 



(31) 



rn 



For Cen A we use a proton jet luminosity of L p = 3 x 10 44 erg s _1 , E m i n = 10 GeV, E max = 10 7 GeV, as determined 
from the interpretation of the Fcrmi-LAT observations reported in Rcf. J]. The L p value is only 10 times larger than 
the kinetic power we used for electrons. If we compare &2 to k e we see that the former is down by a factor of m^/m^ 
with respect to the latter. This means that for protons we loose about six orders of magnitude in the normalization 
of the energy distribution, compared to electrons. We thus generically expect a much lower signal here than in the 
electron case. Recall that with the electrons the kinematic condition on the minimum energy of the incoming electron 
read: E m i n = — E — 3 . The same condition applies to the incoming quark here, and it translates into a condition 



on x m in, the lower limit of the integral in Eq. (27) 



x . - E il E *> f 32 ) 

l-^(l-COS0) 

Also, we must have x m i n < 1, that sets an upper limit on the photon energy 

El < l/E p + l/M x (l-cos6) (33) 

Similarly to the electron case, when the minimum energy of the incoming quark is such that we lose the first resonance 
(s > M|), we expect the signal to drop. This happens for 



M y (M? - Ml) 
7 2M2 + (Mf - M»)(l - cos ' 



We show in Fig 12 our results for proton jets in the Cen A AGN, which confirm the theoretical expectations 
outlined above. In the Figure, we set 5dm = 10 s M /pc 2 and L p = 3 x 10 44 erg/s -1 . The green, blue and purple 
lines correspond to mono-energetic proton jets with increasing energies set to E p — 10 2 , 10 3 and 10 4 GeV, while the 
brown line refers to an E p proton energy spectrum. Note that when the masses of neutralino and of the squarks 
(which, similarly to the electron, we here assume to be all degenerate in mass) are close to each other (as in the right 
panel), the signal drops by a relatively small amount, as it is still enhanced by the second resonance, s' = M|, that 
depends on w, until it falls down again at even larger energies because the integral over x is restricted to regions of 
high x, where the PDF drops. In the plot on the left, for E p = 100 GeV we never hit the first resonance. This is the 
reason why the signal is lower than for all other assumptions on the proton spectrum, and why the signal drops down 
smoothly, simply due to the proton PDF. 

The Figures suggest that in the case of proton-dominated jets, the detection of photons out of proton-dark matter 
scattering is always several orders of magnitude smaller than the gamma-ray output from nearby AGN jets as detected 
by Fermi-LAT and by H.E.S.S. The latitude one has in the overall normalization factors does not seem to allow for 
an enhancement as large as it would be needed to get a detectable signal. While in Fig. [12] we consider only the case 
of Cen A, we also verified that the situation is even less promising for M87. 



VII. CONCLUSIONS 



Recent data from the Fermi Large Area Telescope and from H.E.S.S. on the flux of gamma rays from the local AGN 
in Centaurus A and M87 motivate the present critical re-assessment of the possibility of detecting the scattering of 
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FIG. 12: The spectral energy distribution for photons resulting from proton-dark matter scattering in the Centaurus A AGN jet, 
under the assumption that the latter is dominantly composed of protons. In the Figure, the dark matter density normalization 
factor 5dm = 10 s Mq/pc 2 and the proton jet luminosity L p = 3 x 10 44 erg/s -1 . The green, blue and purple lines correspond to 
mono-energetic proton jets with increasing energies set to E p — 10 2 , 10 3 and 10 4 GeV, while the brown line refers to an E 2 
proton energy spectrum. The Fermi-LAT data are from Ref. [5], while the H.E.S.S. data from Ref. [3]. 



high-energy electrons in AGN jets off of dark matter. This novel indirect dark matter detection method was originally 
proposed in Ref. £Q where, however, for several reasons reviewed here, pessimistic conclusions were reached. In the 
present study we re-evaluated both the dark matter density distribution around compact objects in the light of the 
new theoretical results of [TT] and [H] and of the observational results of [HI EI]; we modeled the AGN jets after the 
significant new observational results of Ref. jH [5] ; finally, and more crucially, we carried out the complete calculation 
of the relevant electron (and proton) scattering cross section off of dark matter. 

We presented the full calculation for the case of supersymmetric neutralinos, and we presented an approximate, 
but, we argued, very accurate, result for the lightest Kaluza-Klein particle of Universal Extra Dimensions. The most 
important result that this calculation produced was the realization that the scattering process can be resonantly 
enhanced. Not only does this resonance produce much larger scattering rates than in the approximate dimensional 
analysis of [T], but it also dictates a very distinctive spectral feature, corresponding to the kinematic limit where the 
scattering can no longer occur resonantly. The spectral feature manifests in the final photon spectrum as a dramatic 
cutoff to an otherwise very hard spectrum. 

We argued that in view of the novel results of the present re-assessment of the dark matter search strategy outlined 
in PP, a large enough signal might be expected, especially for suitably light dark matter particle masses and for low 
enough masses associated to the dark sector partner of the electron (i.e. the selectron for the supersymmetric case 
and the Kaluza-Klein electron for UED). We compared our predictions with the actual gamma-ray data collected by 
Fermi-LAT, and found that for reasonable assumptions on the AGN jet energy content, spectrum and composition, a 
signal from electron-dark matter scattering at the level of the detected photon flux can be expected. We showed that 
photon self-absorption is irrelevant for the attenuation of the resulting gamma-ray flux, and that if the AGN jets are 
proton- instead of electron-dominated, the expected photon flux is likely much below the detectable intensity. 

While the complicated nature of the gamma-ray emission from AGNs naturally prevents any conclusion on the 
possibility that the detected photons be exotic and possibly resulting from the interactions of high-energy electrons 
in the jet with the dense dark matter environment, it is exciting that such a signal might in principle be detectable. 
For Centaurus A, in particular, the difficulty, pointed out in Ref. [5], in explaining the low energy and the high 
energy emissions with a single-zone synchrotron/synchrotron self Compton (SSC) model leaves some room to explore 
other possibilities, such as the one proposed in this study. In the optimistic scenario that future collider results point 
towards a dark sector spectrum with features similar to those expected to produce a signal in AGN jets, further 
studies of the process examined here will most certainly deserve considerable attention. 
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Appendix A: Particle energy distributions in the AGN jet boosted frame 

We consider a model for the AGN jet consisting of a blob moving along the z axis in the black hole rest frame with 
a bulk Lorentz factor Y b = (1 — Pb)^ 1 ^ 2 ■ The axis of the blob frame is assumed to lie along the same direction 
as z, 9' and 9 are the polar angles in the blob and in the black hole frame, respectively. We define = cos 9' and 
fx = cos 9. Primed quantities refer to the blob frame, unprimed quantities to the black hole frame. 

Consider electrons moving isotropically in the blob rest frame with the following power law distribution 



W$AGN | / / 

H^(y)d7V = ifceY-* i+ T 

d l 2 V7br 



dj'dn'. (Al) 



Our goal is to boost this spectral distribution and calculate it as resulting in the black hole frame. The components 
of the electron velocity parallel and perpendicular to the z axis transform, respectively, as 

= [3 cos 9 - p B Rl = ftsinfl 

P|1 l-/3 B /3cos0' p± r B (i- pBpcosey K ' 

and it is straightforward to compute 

V = (1 - P\\ - P'lY 1 ' 2 = (1 - PbP^Tb. (A3) 
Since the electrons are highly relativistic, we can safely approximate /3 ~ 1 in the above and write 

y = (1 - Pb»)iTb. (A4) 
With this approximation, the polar angle transforms as 

M' = f^- (A5) 
1 - PbM 



Working out the differentials for the two Equations above is easy, and we can then simply substitute in Eq. ( Al I to 
find 

j $ agn d$ AGN du 

—^y-h'W d ^ = -^hrsii - pBtJ)) (T B (1 - Mdri) ^ _^ . (A6) 

It is instructive to look at the plot of /i' as a function of ii, which we show in Fig. [Ajfor Tb = 3, to appreciate how 
most of the electrons are seen to move close to the forward direction. 



Appendix B: The differential cross section for the supersymmetric case 

,2 

In this appendix we compute the differential cross section a dfl for the process "electron + neutralino — > electron 
+ neutralino + photon" in the context of the supersymmetric extension of the Standard Model, and for degenerate 
selectron masses. As discussed above, when the exchanged selectron is on-shell, the only relevant Feynman diagrams 
are the ones shown in Fig. [2j 

1. Bremsstrahlung form the initial electron: 

Ml = telu s/ e'(a L P R + a R P L ) X 'x(¥-^(aLPL + a R P R )e, (Bl) 
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FIG. 13: The polar angle [i = cos#' in the blob frame versus the corresponding angle fj, — cos 6 in the black hole frame. 



where (a R ) is the coupling for the vertex left (right)-handed electron - selectron - neutralino 4 , and Pl,r = 
|(1 -F 75) are the corresponding projectors. Also, the identity (e£) = eP R was used. 

2. Bremsstrahlung form the final electron: 

M>i = ie^U s e , {a L P R + a R P L )^(l/ + q^x' x(a L P L + a R P R )e. (B2) 

3. Photon emission from the selectron: 
M% = ieU s U s ,(p + k+p' + k'Y (l-i T — ) e'{a L P R + a R P L ) X 'x{aLPL + a R P R )e. (B3) 



It can be shown that the full amplitude is gauge invariant, ^2 i=1 2 3 Qn-^-i = 0- To this end, one needs the definitions 

given above, and the identity q a {p + k)^ — q^ip' + k') 1 * = Note that the Feynman rule in the last amplitude 

is a generalized one for unstable particles, that was introduced to make the amplitude explicitely gauge invariant. It 
reduces to the usual one for T — 0. To calculate the spin-averaged squared amplitude we first consider terms with a 
1/t' divergence, which will give the leading logarithmic terms. These are 



\ £ \e;(\)MZ\ 2 = 2e\ai + a R )^(pkWq), 

A, spins 

\ J2 K(A)^(A)^r+w = - 2e2 (4+4)^in S 'i 2 [s / -Mi + ^ r 

\, spins 

\ {e;(\)M%e v (\)Mr + h.c} 



}(p'k')4(pk) 2 , 



4 

A, spins 



2e 2 (al + af^Re(IIJI 



(pk)(p'k') + (p'k')ipk') + [pk){kk') - { -^[A{p'k'){pk) + (pq)(s - s')} 



, (B4) 



where we used t' = (k' + q) 2 = 2(k'q) and we neglected (k'q) in the numerator where possible. For future convenience 
we can group them into 

\ M \L=l E KW^2\ 2 + \ E KWM^(\)Mr + h.c} + 1 - J2 {e;(\)M%e v {\)M? + h.c}. (B5) 

A, spins A, spins A, spins 



4 a? L = 2g 2 tan 2 8\y and a? R = \a? L , with g the Standard Model SU(2) gauge coupling and Oyy the Weinberg angle. For the derivation of 
these couplings in the MSSM see, e.g., the appendix of Ref. 1241 . 
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The remaining pieces that do not have large logarithmic enhancements are 
\ ^ KWM?f = 2e 2 (ai + aj i )^(p'k')( Pq ), 

A, spins 

\ X |<(A)A^| 2 = -2e>{a\+a%)\Il s \ 2 \n sl \ 2 {pk){p'k'){s + s') 

\,spins 



s + vV) 



- £ {e;(X)M^(X)Mr + h.c} = -2e\a\ + 4)M|II s ,| 2 [.s - M? + ]{p'k')[A{pk){p'k') - g] 



A, spins 



Again, we group them into the following term: 



A, spins 

The differential cross section is then 

\4r4 



X, spins 



A, spins 



da = (27r) 4 <5 4 (p + fc - p' - k' - q) - 



1 



d 3 k' 



d 3 q d 3 p' 



AM X E 2£'(2tt) 3 2^ 7 (27r) 3 2E' n (2tt) 
The kinematics were discussed in section [Tvl Note that, because the scalar products 



(\M\l R + \M\l olo , 



(B6) 



(B7) 



(B8) 



{p'k') = (p+k-k' -q,k')^ {pk) -(p + k, q) 
(p'q) = (p + k - k' - q,q) -> (p + fc, g) 



(B9) 



do not depend on fc', the integral over the solid k' angle is trivial for most terms in the squared amplitude. The final 
result for the cross section in the supersymmetric case is 



f 



dE 7 dQ (2tt) 5 32M x ^ n 
tt 1 



X —r (\M\f os t'J dn k ,^+4n\M\l olo ^j 



(2tt) 5 32M x EE n 



r [\M\l g t'ln 



4E n 



4\M\ 2 



nolog 
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